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In the field theoretical description of hadronic scattering processes, single transverse-spin asym- 
metries arise due to gluon initial and final state interactions. These interactions lead to process 
dependent Wilson lines in the operator definitions of transverse momentum dependent parton dis- 
tribution functions. In particular for hadron-hadron scattering processes with hadronic final states 
this has important ramifications for possible factorization formulas in terms of (non)universal TMD 
parton distribution functions. In this paper we will systematically separate the universality-breaking 
parts of the TMD parton correlators from the universal T-even and T-odd parts. This might play 
an important role in future factorization studies for these processes. We also show that such factor- 
ization theorems will (amongst others) involve the gluonic pole cross sections, which have previously 
been shown to describe the hard partonic scattering in weighted spin asymmetries. 



I. INTRODUCTION 



Many theoretical as well as experimental studies in recent years have been aimed at better understanding the 
processes that cause spin asymmetries in hadronic scattering. A mechanism to generate single-spin asymmetries 
(SSA) through soft gluon interactions between the target remnants and the initial and final state partons was first 
proposed in the context of coUinear factorization P, H, H, 0, S B 0, 11] ■ This coUinear factorization formalism involves, 
apart from the usual twist-two quark correlators, also twist-three coUinear quark-gluon matrix elements. Since they 
contain the field operator of a zero-momentum gluon, they are referred to as gluonic pole matrix elements. An 
important example is the Qiu-Sterman matrix element Tf(x,x) [J, 0, S 0, 0] ■ 

Several other mechanisms to generate SSA's through the effects of the intrinsic transverse momenta of the partons 
have also been proposed. For instance, in the Sivers effect the asymmetry arises in the initial state due to a correlation 
between the intrinsic transverse motion of an unpolarized quark and the transverse spin of its parent hadron [loj . 
The effect can be described by a transverse momentum dependent (TMD) distribution function f^{x,p'^) pd| . Such a 
function can exist by the grace of soft gluon interactions between the target remnants and the active partons [13, [3] . 
These interactions give rise to process dependent Wilson lines, also called gauge links, in the definitions of TMD 
parton distribution and fragmentation functions. The Wilson lines secure the gauge invariance of these definitions. 
At the same time they prevent the use of time-reversal to argue that the Sivers effect should vanish. Instead, time- 
reversal can be used to derive non-trivial 'universality' relations between the Sivers functions in different processes. 
For instance^Jt was shown that the Sivers function in SIDIS, which contains a future pointing Wilson line, has opposite 
sign [H, as the TMD function in Drell-Yan scattering, which involves a past pointing Wilson line. Moreover, 

the Wilson lines are also crucial ingredients in the derivation [T5| of the relation between the Sivers function and the 

Qiu-Sterman matrix clement, 2M f^^^\x) = —gTp{x,x), demonstrating that the first transverse moment of the Sivers 
function is a gluonic pole matrix element. 

The process dependence of the Wilson lines in TMD parton correlators makes the study of the (non)universality 
of these functions particularly important. In the basic electroweak processes, SIDIS, Drell-Yan scattering and e+e"- 
annihilation, the hard partonic parts of the process are just simple electroweak vertices (at tree- level). Depending 
on the particular process only initial or final-state gluon interactions contribute and, correspondingly, only future 
and past pointing Wilson lines occur. However, when going to hadronic processes that involve hard parts with more 
colored external legs, such as in hadronic dijet or photon-jet production, there can be both initial and final state 
gluon interactions. As a result, the Wilson lines resulting from a resummation of all exchanged coUinear gluons will 
also be more complicated than just the simple future and past pointing Wilson lines [l6l. [TtI Ii8|. In particular, for 
each of the Feynman diagrams that contribute to the hard partonic part of the hadronic scattering process there is, 
in principle, a different gauge link structure. 
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For the TMD distribution functions this at first sight seems to comphcate things considerably. However, for the 
colhnear distribution functions remarkable simplifications occur. Upon integration over intrinsic transverse momenta 
all the effects of the complicated gauge link structures in the TMD correlators disappear, while for the transverse 
moment they contribute a gluonic pole matrix element with multiplicative prefactors, referred to as gluonic pole 
strengths. These are color-fractions that, in principle, differ for each Fcynman diagram that contributes to the 
partonic subprocess. Therefore, for a given subprocess one can multiply the color factors with the contribution of 
each partonic diagram and collect them in modified (but manifestly gauge invariant) hard cross sections fvA [l^ . 
These modified hard functions, called gluonic pole cross sections, appear whenever gluonic pole matrix elements (such 
as the first moments of the Sivers and Boer-Mulders functions) contribute. This is typically the case in weighted 
azimuthal spin asymmetries. 

The effects of the gluon initial and final-state interactions for the fully TMD treatment of these processes is less 
clear-cut. In Refs [lo, [2l| a TMD factorization formula based on one-gluon radiation was proposed for the quark- 
Sivers contribution to the SSA in dijet production in proton-proton scattering. This result involves the gluonic 
pole cross sections found in Refs [13, [Si hard parts, folded with the TMD distribution functions as measured in 
SIDIS {i.e. with a future pointing Wilson line in their definitions). On the other hand, in Refs [IB, [13, HBl it was 
observed that complicated Wilson line structures occur in the TMD distribution (and fragmentation) functions in 
such processes. Those results, in concurrence with a model calculation, led the authors of Ref. [l^ to conclude that 
a TMD factorization formula for spin asymmetries in processes such as dijet production in proton-proton scattering 
cannot be written down with universal distribution functions. It is also asserted that a proof of TMD factorization 
for such processes will be essentially different from the existing proofs for SIDIS and Drell-Yan scattering and that 
it will probably involve 'effective' TMD parton distribution functions [1^. Recent extensions (23, '25l| of the work 
in [13, Hlj, [13] also include the contributions of two colhnear gluons (as was previously discussed for Drell-Yan [26j). 
These indicate that the Feynm an g raph calculations in Refs [la, [13, E, [l3| , [13, HH and [l^l are "mutually consistent" 
up to two-gluon contributions [24 1 . 

By using the gluonic pole strengths we will in this paper systematically separate the universality-breaking parts of 
the TMD parton correlators from the universal T-even and T-odd matrix elements. It is a non-trivial observation 
that this is possible and we believe that it constitutes another important ingredient in trying to relate the results 
of Refs [H, [13, [H, [i3| and Refs [13, 113] . We demonstrate that the gluonic pole cross sections are also encountered 
in unintegrated, unweighted processes. In particular, we will argue that the gluonic pole cross sections can also 
emerge in unweighted spin-averaged processes and that ordinary partonic cross sections can also arise in unweighted 
single-spin asymmetries, though they appear in such a way that they will vanish for the integrated and weighted 
processes [13, [l3, [13, [Hi, respectively. We will start by recapitulating the colhnear case in sectionlllland the appearance 
of universal colhnear correlators in hadronic cross sections in section IIIII The study of the non-universality of the 
TMD parton correlators will be presented in sections IIVI and |Vl followed by a discussion on how the non-universal 
TMD correlators affect hadronic cross sections (section |VI|. After summarizing in section [VTIl we list all universality- 
breaking matrix elements that are encountered at tree- level in 2^2 hadronic scattering processes (appendix [X)) . 



II. COLLINEAR CORRELATORS 



For a twist analysis of hadronic variables in high-energy physics it is useful to make a Sudakov decomposition 
pf^ ~ xP^+an^^+pl^ of the momentum of each active parton. The Sudakov vector n is an arbitrary light-like 
four- vector n^ = Q that has non-zero overlap P-n with the hadron's momentum P^. We will choose the Sudakov 
vector such that this overlap is positive and of the order of the hard scale. Up to subleading twist its coefficient 
a — {p-P—xM'^)/ (P-n) is always integrated over. The vector pT is called the intrinsic transverse momentum of the 
parton. It is orthogonal to both P and n, i.e. pt-P =PT'n = 0, and will appear suppressed by one power of the 
hard scale with respect to the colhnear term. Vectors in the transverse plane can be obtained by using the transverse 
projectors = g^^^ —P^'^n'^^ / (P-n) and e^" = t'^'^P'^ PpU^r / {P -n) . Note that each observed hadron can have a different 
transverse plane. 

We consider colhnear quark distribution functions as being obtained from transverse momentum dependent (TMD) 
quark distribution hmctions. Those are projections of the TMD quark correlator defined on the light-front (LF: 
C-n = 0) 

^f[x,pT;P,S) ^ I '^^^'^2P^ "^"'^ (P,5|^,.(0)^o;«] MO |P,5)Jlf • (1) 

The Wilson line or gauge link ^,,;^] = 'Pexp[— igj^ ds-A"(s) f° ] is a path-ordered exponential along the integration 
path C with endpoints at rj and ^. Its presence in the hadronic matrix element is required by gauge-invariance. In the 
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FIG. 1: Simplest structures (without loops) for gauge links and operators in quark correlators (a)-(b) and gluon 
correlators (c)-(f). 



TMD correlator ([T]) the integration path C in the gauge link is process-dependent. In the diagrammatic approach the 
Wilson lines arise by resumming all gluon interactions between the soft and the hard partonic parts of the hadronic 
process [l^ [H, [29l Isol . Consequently, the integration path C is fixed by the (color-flow structure of) the hard 
partonic scattering [18"|. Basic examples are semi-inclusive deep-inelastic scattering (SIDIS) where the resummation 
of all final-state interactions leads to the future pointing Wilson line U^^^ , and Drell-Yan scattering where the initial- 
state interactions lead to the past pointing Wilson line see Figures [TJi and b. All Wilson lines in this text are in 
the three-dimensional fundamental representation of the color matrices. 

Going beyond the simplest electroweak processes such as SIDIS, Drell-Yan scattering and e^e^-annihilation, the 
competing effects of the gluonic initial and final-state interactions lead to gauge link structures that can be quite 
more complicated than the future or past pointing Wilson lines [l6l . [itI [18| . The situation becomes particularly 
notorious when considering processes which have several Feynman diagrams that contribute to the partonic scat- 
tering. In that case each cut Feynman diagram D can, in principle, lead to a different gauge-invariant correlator 
$["l(a;,pT) = 4''"^^^l(a;,PT) = ^''^l(2;,PT) [Hi- This observation leads to a broad spectrum of different TMD 

parton correlators that appear in hadronic scattering processes. 

For coUinear correlators the situation is simpler. For instance, in the pT-integrated correlator defined on the 
lightconc (LC: ^■n = ^T = 0) 
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all process-dependence of the gauge link disappears, leaving just a straight Wilson line U^.^^ in the lightconc n- 
direction, where n is the lightlike vector in the Sudakov decomposition of the quark momentum p (we will use the 
non-calligraphed letter U to indicate straight line segments). Another situation is encountered in the transverse 
momentum weighted correlators (the transverse moments). In the transverse moments a (sub) process-dependence 
remains as a direct consequence of the presence of the gauge links in the TMD correlators. Nevertheless, a simple 
decomposition can still be made (omitting the Dirac indices) (isl . [TtI . [Tgj : 
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The only process dependence due to the Wilson lines in the TMD correlators resides in the multiplicative factors 
= Cq^^^^ = Cq'^ . They are color-fractions that are fixed by the color-flow structure of the hard partonic function 
of the scattering process [l7l.[l9l|. We will refer to them as gluonic pole strengths. Important examples are the transverse 

moments of the correlators $^+1 in SIDIS and <i>t^l in Drell-Yan scattering, for which one has C^^' = Cq ' = ±1 p^ . 
Transverse momentum dependent gluon distribution functions are projections of the TMD correlator 

T^^^^'^'^-^{x,pt;P,S) = ^ I ^^^'g^f^^ e^^-«Tr(F,5|F''''(0)W[o^,]F-(0%o] I^,^)Jlf ■ (6) 

Here Tr indicates a trace over color-triplet indices. Writing the field-operators in the color-triplet representation 
requires the inclusion of two Wilson lines ^o;^] ^-nd ^j'^.q] [ist - They again arise from the resummation of gluon initial 
and final-state interactions. In general this will lead to two unrelated Wilson lines U and W . In the particular case 
that W —U\ the gluon correlator can also be written as the product {FaUabFb) of two gluon fields with the Wilson 
line U in the adjoint representation of SU{N). This is for instance the case for the gluon correlators in Figs[T]; and d, 
but not for the gluon correlators in Figs [TJj and f . 

In the pT-integrated correlator on the lightcone the process dependence of the TMD gluon correlator disappears, 

r'^-{x) = Jd^PT r['^^"'l^''(x,pT) = ^^|^^e-«-^)Tr(F,5|F^''(0)C/f^^^]F^^ (7) 

However, as for the quark correlator, a subprocess-dependence due to the Wilson lines in the TMD gluon correlators 
remains in the transverse moments. The analogue of the decomposition ([3]) in the case of the gluon correlator is (with 
r["'"'l(a;,]3T) =r["(^)'"'(^)l(a;,pT) = r[^l(a;,pT) and omitting the gluon field indices n and i^) 

r[fi"(x) = nix) + ^v^G.M + 4" ^rg^(x,x) . (8) 

The matrix elements ^Gf and Tq^ are the two gluonic pole matrix elements that correspond to the two possible 
ways to construct color-singlets from three gluon fields [l^, IMI ■ They involve the antisymmetric / and symmetric d 
structure constants of SU (3) , respectively. The only process dependence coming from the Wilson lines in the TMD 
correlators is contained in the gluonic pole strengths Cq^'^^ = C^^'^^ \u{D)m {d)\ ^ rpj^^ collinear correlators are 

ix' (ri-P) i(x — x'){^-P) 
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and 



T%{x) = Tf,{x) - / dx' P- V^aM^x-^') ■ (10) 



The collinear (anti)quark and gluon fragmentation correlators can be analyzed in the same way. The matrix elements 
in ([3]) and ((HI) contain the collinear T-even and T-odd parton distribution functions, see e.g. [19|. 



III. COLLINEAR FUNCTIONS IN HADRONIC CROSS SECTIONS 



In the diagrammatic approach, the calculation of the hadronic cross sections starts off with the transverse momentum 
dependent parton correlators (TMD distribution and fragmentation functions), which will appear in combination with 
squares of hard partonic amplitudes. In general the hard amplitude contains more terms, that is H = Hi. In the 
squared amplitude one therefore has terms like F[*Hj = E^-^l, where D refers to the cut Feynman diagram that is the 
pictorial representation of the product of the amplitude Hj and conjugate hard amplitude H* . The hadronic cross 
section da of a hadronic scattering process mediated by a two-to-two partonic subprocess a{pi)h{p2)^c{ki)d{k2) and 
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where the outgoing hadrons and/or jets are produced with large perpendicular component with respect to the beam 
will contain the following structure in the integrand: 



,d D 



where the parton momenta are approximately (compared to the hard scale) on-shell. The convolutions '(8)' represent 
the appropriate Dirac and color traces for the hard function Sl^l. To get to the hadronic cross section one has 
to multiply by the flux factor and integrate over the final-state phase-space and parton momenta including a delta 
function for momentum conservation on the partonic level. Since our aim in this paper is to display some general 
features that a /cr-factorization formula (if it exists) will have due to the process-dependence of the Wilson lines that 
arise in the diagrammatic TMD gauge link approach, we focus our discussion on the Wilson lines and neglect soft 
factors. In the full fcr-factorization formula such factors will most likely also be present to account for soft-gluon 
effects. 

The TMD correlators in (jlip are the gauge invariant non-universal (anti)quark/gluon correlators that contain the 
appropriate Wilson lines for the particular color-flow diagram D that contributes to the partonic subprocess ab^cd. 
This is the reason why in expression pip the summation over cut diagrams D is displayed explicitly. The hard 
functions, i.e. the expressions S^^l of the individual Feynman diagrams are, themselves, not gauge invariant. For 
azimuthal dependence originating from only one of the partons one can effectively use the correlators calculated in 
Ref. Furthermore, in the tree- level discussion employed here the Wilson lines are along the lightlike n-direction, 
though a non- lightlike ^0 direction may be required when higher-order corrections are taken into account [33 , [33| . 

From momentum conservation on the partonic level it will follow that, depending on the process, some components of 
the partonic momenta can be measured {e.g. in a way similar to the identification of the incoming parton momentum 
fraction x with the Bjorken scaling variable xb in deep inelastic scattering). This also works for the transverse 
momenta. For instance, for a hadronic scattering process with a two-to-two hard subprocess the structure in (jlip will 
appear with a delta function for momentum conservation enforcing the relation pi+P2^ki—k2 = 0. There are ways to 
measure one or several components of qr =piT+P2T—kiT—k2T ^ Ki/ Z1+K2/ Z2—X1P1—X2P2, which is not required 
to vanish by momentum conservation since the directions of the intrinsic transverse momenta of the partons can be 
different for each observed hadron (in back-to-back jet production in hadron-hadron scattering it is the component 
along the outgoing jet direction in the plane perpendicular to the beam axis that is experimentally accessible through 
the relation q^-K^^ oc sm{S<j)), where S<j) is the azimuthal imbalance of the two jets in the perpendicular plane [TtI. uM. 
[13, [13]). This quantity defines a scale much smaller than the hard scale of the process. Using these components one 
can construct integrated and weighted hadronic cross sections. Integrated cross sections will involve the structure 

^{xi,X2,Zi,Z2) = Jd^PiTd^P2T(fkiTd^k2T T,{pi,p2,ki,k2) 

= E Tr { $a(xi)$6(x2) tab^cdixi,X2,Zi,Z2) Ac(zi) Ad(z2) } , (12) 

while weighted cross sections will involve (making use of the decomposition in Eq. ([3])) 

T,^q{xi,X2,zi,Z2) = j d^p^j.d^P2Td^k^j.d'^k2T Pit 2(^1,^2,^1,^2) 

= E Y.^^{^^ad'"i^i)'^b{x2) S[^l(a;i,a;2,zi,Z2) A,(zi)A,(z2) } 

a,---,d D 

= E [ Tr{ $;ja(a;i)$h(a;2) Sa6^cd(a;i,a;2,zi,Z2) Ac(zi)Arf(z2) } 

a,- ■ ■ 

+ Tr { 7r$"Q(a:i,xi)$fc(a;2) E[a]fc^cd(2^i,a;2, 2:1,22) Ac(zi)A<i(z2) } ] , (13) 

and similar expressions £20, '^v d and S2'a which are obtained by weighting with P2T, fciT and k2T, respectively. In 
these expressions only universal coUinear correlators appear. In Eqs and pH]) we have defined the hard functions 



Sa6^cd(a;i,a;2,zi,Z2) = i;[-°l(a;i,a;2, 2:1,22) , (14a) 



^[a\b^cd{xi,X2,Zi,Z2) {xi,X2,Zi,Z2) . (14b) 
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The factors Cq (a) are the gluonic pole strengths that appear in the decomposition of the transverse moment of the 
TMD correlator of parton a. In expression (|13p this parton was implicitly taken to be a quark. If it were a gluon there 
would have been two S[g]f,^cd terms, one corresponding to each of the gluonic pole matrix elements and Tq^, cf. 
Eq. ^ or Ref. [l^. The hard functions in Eqs (|14ap and (|14b[) no longer depend on the individual (diagrammatic) 
contributions D, but only on the hard process ab^cd. Moreover, in contrast to the hard functions S^^l that appear 
in (jlip . they are gauge invariant expressions. After performing the traces the £ab->cd reduce to the partonic cross 
sections daab^cd and the S[a]b^cd reduce to the gluonic pole cross sections d(7[a]b~>cd calculated in Refs [itI. ITqI. [28l| . 



IV. TRANSVERSE MOMENTUM DEPENDENT CORRELATORS 



To study aziniuthal asymmetries arising from one of the partons in hadronic processes mediated by 2— >2 partonic 
subprocesses at tree-level it is possible, as we will show explicitly, to organize the TMD correlators in a decomposition 
analogous to ([3]) containing TMD correlators with special properties: 

$[^l(a;,pT) = <P^'"'-'''\x,pt) + C^f 1 ^$(,"'^^''(x,pt) . (15) 

Here D refers to a particular cut Feynman diagram that contributes to the cross section of the partonic process 
ab^cd. The gluonic pole factors C^' are the same as those in the decomposition of the collincar correlators in ([3]). 
An important difference between the decomposition of the coUinear correlator in ([3]) and the decomposition of the 
TMD correlator in ()15|) is that the matrix elements in the latter decomposition are not universal in general. They 
depend on the partonic process ab^cd but, in contrast to the TMD quark correlators <I>t^l on the l.h.s. of the 
decomposition, they do not depend on the individual cut Feynman diagram D. The only diagram dependence resides 
in the gluonic pole factors. The matrix elements on the r.h.s. of ([T5|) have been chosen such that they reduce to 
the familiar universal (process independent) coUinear matrix elements when integrating over or weighting with the 
intrinsic transverse momenta: 

d'^PT ^^''^~""^Hx,pt) = <^{x) , J d'^pj, p^ ^^''''-""^\x,pt) = ^oix) . (16a) 

J d'pT ^^f-'''\x,pT) = , I d'p^ p^ ^^f-^'\x,pT) = ^U^,x) . (16b) 

The most straightforward illustration of the decomposition (fT5|) for quark TMDs are the quark correlators {x,pt) 
in SIDIS and ^^^\x,pt) in Drell-Yan scattering, which contain the simple future and past-pointing Wilson lines W^^l 
and ^[^1 , respectively. For those correlators one has [15| 

$[±l(a;,pT) = $(^-'=™")(x,pt) + Cl^^ $(^-°'^^)(x,pt) , (17) 
with the T-even and T-odd quark correlators 

$(^-°™") {x,pt)^^{ $ W {x,pt) + $[-1 {x,pt) ) , (18a) 
{x,pt) = \ ( $ W (x,pt) - $[-1 {x,pt) ) . (18b) 

The factors C^^' =±1 are the same as those for the transverse moments ([3]) in those processes. 

In contrast to p7|) . we observe that the TMD matrix elements ^'^°'^^'^'^\x,pt) and 'k^'^q'^'^'^^ {x,pt) on the r.h.s. 
of in general do not have definite behavior under time-reversal. However, the process dependent universality- 
breaking parts of the TMD correlators can be separated from the universal T-cven and T-odd parts in (|18ap and 
(Il8b|: 

i:<^f~''^\x,PT) = $(^-°'^^)(x,pt) + ^8^f-''\x,pT) . (19b) 

In these expressions all process dependence due to Wilson lines on the light- front is now contained in process-dependent 
universality-breaking matrix elements (5<I>(°^^'^'')(x,pt) and 7r(5$g''^^'^'' (a;,pT), which we will refer to as junk-TMD. 
Also these in general have no definite behavior under time-reversal, but they do have the special properties that they 



7 



vanish after pT-integration and weighting: 

J dPpT 6^^''''^'"^\x,pt) = J d?PT Pt 5^''''^^"^\x,pt) - , (20a) 
d^T S^f-'^Hx^pT) - / d^pr p^ S^'^f-''\x,pT) = . (20b) 



i-G \-^iFT) — j FT PT "^G 

This is consistent with the properties expressed in p6)) . The expressions in Eq. (|20p will be used in section IVII 
to show that the universality- breaking matrix elements vanish in integrated and weighted hadronic cross sections. 
Moreover, as can be seen from their explicit expressions in appendix [Al the (anti)quark/gluon universality-breaking 
matrix elements vanish in an order g expansion of the Wilson lines {i.e. the one-gluon radiation contribution). 

All universality-breaking matrix elements that occur at tree-level in proton-proton scattering with hadronic final 
states are listed in appendix]^ The TMD correlators '^^^\x,pt) in these processes have already been derived in 
Ref. [HI and are given in the tables of that reference (the TMD correlators and gluonic pole factors that appear at 
tree-level in direct photon-jet production in proton-proton scattering can be found in Ref. [35]). It is straightforward 
to verify that these results are reproduced with the matrix elements given in appendix [X] and through the decompo- 
sitions ()15p and (|19p . This should not come as a surprise, since the matrix elements in the appendix were defined 
that way. It is a remarkable and non-trivial observation that with the gluonic pole strengths all the quark correla- 
tors encountered in a certain partonic process a6— >cd can be decomposed in terms of only the two matrix elements 
^(ab^cd) g^j-^j 7r$[?''^'^'''' (with the properties in Eq. (fT6|) ). It should be mentioned, though, that this decomposition 
is not unique. For instance, one could also have made a decomposition in terms of more matrix elements. That is, by 
including matrix elements that do not contribute to the zeroth (pr-integration) and first (pT-weighting) transverse 
moments in py, but do contribute to the second moment, third moment, etc. It is conceivable that the inclusion of 
these additional matrix elements will allow one to summarize the TMD quark correlators encountered in different 
partonic processes, in the same way as it was possible to summarize all quark correlators associated to the different 
Feynman diagrams D that contribute to one specific partonic process ab^cd by the two matrix elements $(ab^cd) ^^^^ 

TT<^^Q^~*'^'^\ At present this is just speculation, though, and the verification or falsification will require more insight 
into the way that the different Wilson line structures contribute to higher transverse moments. However, regardless 
of all these cautionary remarks we believe that the notational advantage, the points concerning gauge invariance of 
the hard functions that will be addressed in section IVII and the possible role that it could play in relating the gauge 
link formalism to the results of Refs [13, [2l| provide more than enough justification for the decomposition in Eq. (llSp . 
In the case of gluon distributions we start by defining T-even and T-odd gluon correlators (c/. Figs[T}:-f), 



r(^-™)(a:,PT) = \{T^+^+'\x,pt)+T^-'-\x,pt)) , (21a) 

r^^i^"'^ {x,pt) = \ ( pi+^+^i {x^pt) - r["'-'i {x^Pt) ) , (21b) 
Tf^-°'"'\x,PT) = i(r[+'-'i(^,PT)-r[-+^i(x,p^)) , (21c) 

where as in Ref. (TH a gluon correlator is called T-odd if it vanishes when identifying the future and past-pointing 
Wilson lines. Note that in contrast to r(2"-ovon) g^j^j r^Jp^'^'', the correlator r|Jp'^^'' cannot be written as a matrix 
element of two gluon fields with a single Wilson line in the adjoint representation. After pT-integration the T-even 
and T-odd correlators reduce to the universal coUinear gluon matrix elements in the expressions in Eqs ([7]) and ([H]): 

(22a) 
(22b) 
(22c) 

Since there are two distinct ways to construct T-odd gluon correlators, it will follow in the next section that there 
are also two distinct TMD gluon-Sivers distribution functions (c/. Eq. (j28bp ). There is actually also a second way to 
construct a T-even gluon correlator: r''-^"°™"^ = i(rl+'~ 1 -|-r[~'+ 1). However, this correlator is not needed in the 
decomposition ((23|) of the TMD gluon correlators, since the difference between F' (^-even) p(T-ovcn) 

is a matrix 

element that vanishes upon pT-integration and pr-weighting. This difference may, therefore, be absorbed in the 
universality-breaking matrix elements ^r(afc^cd) defined in (|24ap . 

A decomposition resembling the one in ([8]) for TMD gluon correlators can almost be made: 

Fl^l {x,pt) = Ft"^"^^^) (x,pt) + '""l rg^^"^) {x^pt) + do'^ F^?^^^^^' (x^pt) . (23) 



d'^PT 




: r{x) , 


j d^Prp 


Pt 


F(^-'=™")(x,pt) 


= rg(a;) , 


d^PT 




, 


j d^PT 


Pt 


^(T-Odd), ^ 


= vrFg^(a;,a;) 


d^PT 


^(T-Odd) , ^ 


, 


J d^PT 


Pt 


^(T-Odd), ^ 


^nr^^{x,x) 
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This leaves only a specific type of matrix elements with colorless intermediate states unaccounted for (we will return 
to this point in a moment). The TMD matrix elements on the r.h.s. of (|23p only depend on the process and not on 

the particular Feynman diagram D that contributes to that process. The multiplicative factors Cq'^ and Cq^ are 
the gluonic pole factors calculated in Ref . [l^ , the same that also appear in the decomposition ^ of the collinear 

correlator. Under p^-integration and weighting the matrix elements T^"'^^'^'^'' {x,pt) and r[?^^^'^''(a;,pT) have the same 
behavior as T^'^~°^^'^'> [x^px) and t''i^j~°^'^\x ,pt) , respectively. Therefore, one can make the further separation 

p(a6^cd) (-^^p^) ^ p(T-cvc„) _^ ^p(afc-.cd) (^^p^) ^ (24a) 

.T^^'-;'^ {x,pt) = r[ ) {x^pt) + nSr^OrT^ ("'^^) ' (24b) 

in which all process-dependence due to Wilson lines on the light-front has been gathered in the universality-breaking 
matrix elements 6T^"''^^'^'^\x,pt) and 7r(5r|^^^'^'^''(x,pT), which have the special properties that they vanish after a 
PT^-integration and pj^-weighting. 

It is also straightforward to check that through the decompositions in and with the universality-breaking 

matrix elements given in appendix \X\ the TMD gluon correlators in the tables 4, 5 and 8 of Ref. [18] corresponding 
to qg^qg, qg^qg and gg^gg scattering are reproduced. However, in the TMD correlators in the tables 6 and 7 for 
the processes qq — >gg and gg — >qq one does not recover terms of the form {P,S\ Tr [F{^)U^°^ Tr [F(0)Z^Plt] \p^s), 
where — U^^^U^^^'^ . These matrix elements involve colorless intermediate states and they do not contribute to 
the pT-integrated gluon correlators T{x) nor to the first transverse moments t\^\x). They can be included in ([^5]) by 
adding diagram-dependent universality-breaking matrix elements which will not appear in integrated and weighted 
hadronic cross sections. However, it could be that they contribute to the second or higher transverse moments. 



V. PARAMETRIZATIONS OF PARTON CORRELATORS 



At leading twist the parametrizations of the different TMD quark correlators are given by 



+ Sl giL{x,pl) 75f + gMx^p^T) Tof } , (25a) 

'f(^-°^''\x,pT) = i { ^hi{x,p'^) 1^ - fMxA) r } , (25b) 

and 

d<P^^'-^'Hx,PT) - i { Sf[''''-^'\x,p'r)r+kSh[''^'-^''\x,p'r) 75[^T,f] 

+ Sl Sh^l '(x,Pt) 75^771— + <5/iit {x,Pt) 75 



2M M ' ' 2M 

M 



Sl ^g^L^" '{x,Pt) 75f + ^. . ^g^^^" \x,p%) 



i Sh 



±{ab^cd)f 2\ [^T,f] ef^^"^ cr±iab~^cd)f 2 \ m \ /r,c- \ 

(^'^"^-^ ~2M M~^^^ \x^PT)r], (25c) 



with similar parametrizations for the matrix elements ttS^'q^^'^'^^x^pt) containing a (different) set of distribution 
functions jy^I'^l''^'^'')^ etc. The quark distribution functions that appear in the parametrizations (|25ap and (j25b[) 
are the familiar T-even and T-odd (respectively) quark distribution functions as measured in SIDIS. On the other 
hand, the quark distribution functions in (I25cp and in the parametrization of 5^g ^'^^ process dependent. From the 
properties in (|20p one finds that the functions 5hiT and 5giL vanish upon p^-integration, for instance 

(SpT <5/r-^^'^)(x,p|) = 0, (26) 



9 




FIG. 2: Possible behavior of the universal distribution function fi{x,p^) (dashed line) and the 
process-dependent function fi'''^'^'''\x,p^) (solid line) as a function of |Prp. Their difference- 
function, 5f^'^^^'^'^\x,p'^) (dash-dotted line), vanishes upon integration over pr- 



illustrated in Figure [51 Also the functions 6hfj^^\ Shf^^\ Sg[l^, Shf"^^ and Sf^rj;^'' vanish, e.g., 

d'PT ^ 5ftT^^'-^'\x,pl)=Q. (27) 

The same holds for the corresponding functions in the parametrization of tt^^^''^'^'^'' (^.p^)- 

For gluon distribution correlators we use the parameterizations of Ref. ^] with the naming convention discussed 
in Ref. W^: 



+iel^''SL gli^[x,p^) -f ie^" ^^^^ gfri^^Pr) | . (28a) 

1 (■ ptSt PtImc'^} I ■S't{a' ''l 

^ (f/d) y^^PT) - 2x \ M \^^Pt) " ^j^j '^iT K^^Pt) 



-'-^S, htl^^''\x,pl) - htS^f/^\.,p^,) \ . (28b) 

In particular one has two distinct gluon-Sivers distribution functions firf'^^\x,p'^) and fix^'^\x,p'^) correspond- 
ing to the two ways to construct T-odd gluon correlators. Their first transverse moments are the functions 
G^^^''^(i)(a;) = /{^^^^/''^(i)(a;) introduced in Ref. [H. The matrix elements 5r("*^'='')(x,p|,) imdT:5T^(i''^~Jl'"^\x,pl) con- 
tain both T-even and T-odd gluon distribution functions. These functions, however, will vanish under pT^-integration 
or weighting, as for the quark universality-breaking distribution functions (c/. (pS]) and ([211) )■ 

The analysis above can be extended to antiquark distribution correlators in the obvious way. Also the treatment 
of (anti)quark fragmentation correlators is straightforward, with the understanding that in the matrix elements 
^(T-ovon/odd)^^^;,^) the superscripts T -even and T-odd refer to the operator structure in the correlators, each of which 
contain both T-even and T-odd fragmentation functions (this also holds for the gluon fragmentation correlators F). 



VI. TRANSVERSE MOMENTUM DEPENDENT FUNCTIONS IN HADRONIC CROSS SECTIONS 

The advantage of the decomposition of the TMD quark (and gluon) correlators in expression (and ([55)) ) is 
evident: instead of having to list the TMD correlators $[^1 (F^^l) for every cut Feynman diagram of a process, the 
two (or three in the case of gluons) matrix elements on the r.h.s. of those decompositions sufhce for that particular 
process. This observation leads to a reduction in the number of TMD matrix elements that need to be considered. 
From a bookkeeping point of view this makes the treatment of TMD correlators more manageable. What is perhaps 
more important is that the decompositions (jlSp and (j23[) will allow us to show that also unweighted, unintegrated 
hadronic cross sections can be written as products of soft parton correlators and hard partonic functions that are 
separately manifestly gauge invariant, analogous to the coUinear case in expressions and ([T5|) . Moreover, the 
hard partonic functions are the partonic cross sections or the gluonic pole cross sections. This is what will be argued 
in this section. The proper context should be within a transverse momentum dependent factorization theorem. For 
most processes such a theorem does not exist yet. We will therefore take as a starting point the assumption that the 
hadronic cross sections will factorize in a hard partonic function and a soft parton correlator for each observed hadron 
separately, and that gluon initial and final-state interactions will lead to the required Wilson lines. We believe that 
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this assumption is sufficiently generic for our conclusions to be applicable for many hadronic processes, in particular 
to back-to-back dijet or photon-jet production in proton-proton scattering. 

By inserting the decompositions p5|) and ((23|) of the TMD parton correlators into the expression for the unintegrated 
hadronic cross section in Eq. pT|) the parton contribution to a hard 2^2 process becomes 

+ •••}, (29) 

which forms the central result of this paper. Again it is implicitly implied that parton a is an (anti)quark, since 
if it were a gluon there would be two gluonic pole terms rGy^Sj^jj,^^^ and Fg^Ej^j^j^^^^. In (f29|) both terms in the 
decomposition of the TMD correlator of (anti)quark a have been given explicitly, while only the first terms of the 
decompositions (|15p and were used for the other partons. The " contains the other possible combinations 

of the terms in those decompositions and also the contributions where parton a is a gluon. 

The hard functions in expression ([29|) are the partonic and gluonic pole cross sections in ()14ap and (Il4bp (in the 
collinear expansions of the hard functions, which have corrections at order 0{l/s)). Hence, the expression in ([29]) 
demonstrates that by using the decompositions in (fT5|) and ((23| and by introducing the gluonic pole cross sections, 
also the unintegrated, unweighted hadronic cross section can be written as a product of soft TMD parton correlators 
and hard partonic functions that are separately and manifestly gauge invariant. After performing the traces these 
hard functions reduce to the partonic and gluonic pole cross sections calculated in Refs [13, E^j HI] ■ Hence, it follows 
that the gluonic pole cross sections that have been seen [l3, 0, [13, represent the hard partonic functions in 

weighted spin asymmetries already appear in the fully TMD cross sections. This conclusion is consistent with the 
work in Refs [13, [HI, where a TMD factorization formula for the quark-Sivers contribution to single transverse-spin 
asymmetries in dijet production is proposed with the gluonic pole cross sections of Refs [itI [l9j as hard functions. 
However, the work in Refs 20, 21] limits to one-gluon exchange and as a result it obtains the distribution functions 
measured in SIDIS (i.e., transverse momentum dependent distribution functions with a future pointing Wilson line 
in the hadronic matrix elements defining them). In contrast, our expression (|29p involves correlators with process 
dependent Wilson lines in the operator definitions and a factorized form with universal distribution functions would 
only be achieved if, regardless of the appearance of process dependent Wilson lines in their definitions, the TMD matrix 
elements in (j29p are identical for all partonic channels and equal those in semi-inclusive deep inelastic scattering. In the 
present context this translates into a vanishing of all univers ality -breaking matrix elements, for which we see no reason. 
Indeed, universality has recently also been disputed in Ref. [24I, where it is argued that a TMD factorization theorem 
for this process with universal distribution functions is not possible. This is confirmed by an explicit calculation 
including the exchange of two collinear gluons [2^. A recent extension [23| of the work in Refs l20l. [21I inclu ding 
two-gluon exchange also points to non- universality. Moreover, it shows that Refs [itL [l^. |20l. l2l| and |22j | 
are consistent to (at least) that order. We want to emphasize that for a full connection the role of soft factors, which 
have been neglected in the present study, should also be investigated. 

With the parametrizations P5|) and (|28p inserted in P5)) . an expression for the hadronic cross section in terms of 
TMD parton distribution and fragmentation functions is obtained. After performing the traces these hard functions 
reduce to the partonic and gluonic pole cross sections encountered in Refs [13, [la [13, [111 ■ As an illustration we 
consider the contribution to (j29p of an unpolarized quark in an unpolarized hadron, important for unintegrated spin 
averaged hadronic cross sections (summations over parton types are understood): 

dcru ~ { h{x^,p\j.) + 5f[°'^^'"^\x^,p\rr)^ } daab^cd + 5f^^"^''^'"^\x^,p\rj.) da^a\b~,cd ■ (30a) 

Taking as a second example the contribution of an unpolarized quark in a transversally polarized hadron, important 
for unintegrated single-spin asymmetries, one finds 

dar ^ '5/^V^^'''*'(2;i'-Pit) ddab^cd + { ftAxi-,PiT) + ^ftr^"^''^'''^^ (xi^pIt) } da[a]b^cd ■ (30b) 

V ' 

The terms with only universal T-even functions will appear folded with partonic cross sections and the universal 
T-odd functions with gluonic pole cross sections. In addition, there are various universality-breaking functions that 
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appear with partonic cross sections or gluonic pole cross sections. A shorter notation for some of the terms in these 
expressions could have been obtained by not extracting the universal T-even and T-odd parts of the distribution 
functions /i and f^, as indicated by the underbraces in ()30|) . In particular, due to these universality-breaking 
functions the gluonic pole cross sections also appear in the unweighted spin-averaged cross sections pOap and the 
usual partonic cross sections also appear in the unweighted single-spin asymmetries (j30bp . However, in the light of 
the properties in ()26|l - (j27|) it is seen that terms with universality-breaking matrix elements do not contribute to the 
integrated and (/^-weighted (where qt is as defined in section |TTT| hadronic cross sections which are expressed in terms 
of the structures in p2p and : 

integrated: (da) ^ T.{xi,X2) oc fi{xi) d&ab^cd , (31a) 
weighted: (grda) ^ (T^ig + j:2d)ixi,X2) oc f^rl;^\xi) da[a]b^cd ■ (31b) 

For back-to-back jet p roduction in polarized proton-proton scattering {p^ p^jjX) this (in essence) reproduces the 
results of Refs (l7l ll9ll27t. while for photon-jet production {p^ p^jjX) it reproduces the results in Ref. [28j . 

Only if the universality-breaking matrix elements vanish in the unintegrated, unweighted cross sections ([29]) and (j30[) 
does one also in the TMD case arrive at the situation with universal functions only. Otherwise, the non-universality 
of the unweighted processes is important. It will affect the results of experiments that try to look at explicit px- 
dependence or that construct weighted cross sections involving convolutions of TMD functions which upon integration 
do not factorize into transverse moments, e.g. when looking at {sui{(j)hi<t's) ): rather than {P-n-i_ sin{(l)hi(f>s) ) asym- 
metries in SIDIS to extract transversity or Sivers functions. 

VII. SUMMARY 

We have argued that the gluonic pole cross sections, the hard partonic scattering functions that are folded with 
the collinear parton distribution functions in weighted spin-asymmetries, also appear in unintegrated, unweighted 
hadronic cross sections. Assuming as a starting point that the hadronic cross section factorizes at the diagrammatic 
level in a hard partonic function and for each of the observed hadrons a soft parton correlator, we have shown that the 
transverse momentum dependent cross section can be written in terms of soft and hard functions that are separately 
manifestly gauge- invariant. The hard functions are the partonic cross sections or the gluonic pole cross sections. The 
latter do not show up in the integrated cross sections. The soft functions are the TMD parton distribution (and 
fragmentation) functions with Wilson lines on the light-front in their field theoretical operator definitions. These 
Wilson lines can be considered as the collective effect of gluon initial and final state interactions. Since they are 
process dependent, the TMD parton distributions are in general non-universal. By systematically separating the 
universality-breaking parts from the universal T-even and T-odd parts, we arrived at an expression that has soft 
parts multiplying the partonic and gluonic pole cross sections, as was also found in Refs [13, HJ. However, in the 
gauge link approach taken here the gluonic pole cross section can also emerge in TMD spin-averaged processes and 
that ordinary partonic cross sections can also arise in TMD single-spin asymmetries. They appear with universality- 
breaking functions that will vanish for the integrated and weighted processes considered in Refs [13, [13, [13, [13 • They 
also vanish at the level of one-gluon radiation contributions, which corresponds to the order g term of the Wilson 
lines. 

The non-universal terms are well-defined matrix elements. All universality-breaking matrix elements that are 
encountered at tree-level in proton-proton scattering with 2— >2 partonic processes have been calculated and are given 
in the appendix of this paper. In particular, the process-dependent universality-breaking matrix elements S^, etc. 
disappear in the simple electroweak processes with underlying hard parts like qj*—^q and qq^^*. We believe that 
the explicit identification of universality-breaking matrix elements is an important contribution to arrive at a unified 
picture of TMD factorization of hadronic scattering processes. 
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APPENDIX A: UNIVERSALITY-BREAKING MATRIX ELEMENTS 

We list the universality-breaking matrix elements that appear at tree-level in proton-proton scattering. To improve 
readability we employ the schematic notation {2n)-^ J d{^-P)(P^T e^^f (P,S'| V'(0)W l^^'S') Jlf °^ {-^{0)U tp {£_)), 
and similarly for the other parton correlators. For the 2— >2 partonic channels with four colored external legs we only 
list the quark and/or gluon distribution correlators. As can be seen from the tables in Ref. [l^ the TMD correlators 
of the other partons are very similar in structure to the correlators considered here and can straightforwardly be 
obtained by comparing the results below to the tables of the stated reference. In particular the correlators in 
qg^qg scattering are simply obtained by comparing to those in qg—^qg scattering. Similarly, the fragmentation cor- 
relators in qq—fgg (gg—^qq) scattering can be obtained by comparing to the distribution correlators in gg—^qq (qq^gg)- 



97 -^g qq^i 

^i^x,pt) = 7r5$[?^*^«^(a:,pT) = (Al) 
jc^(??-.7*)(^^p^) = tt5^'^^''^^'\x,pt) = (A2) 



S^^^'^-''Hx,pt) = 7v5^1S^-'-'\x,pt) OC (^^(0) { iiy^wWt - iwWt } ^^(0 ) (A4) 
(99-97) (^^^^) ^ 7r(5f [jp"^' {z,kT) = vr^r^^;^'^' {z,kT) = (A5) 



qg^qi 

J$(99->97)(^^p^) = -7r5$(?«-*^^)(x,pT) OC (^(0) { i^^Mw" - i^'^ } ^(0) (A6) 

5^{99->97) ^ tvSA''^/-"'''^ {z,kT) = (A7) 

5r(99-97)(^^p^) OC (Tr [F(0)wWF(o{5W'"'^-5WWt| +7r(0)w[-lF(C){iw'+'^-iw'-'t}] > (A8) 

^5r(?;^'^) {x,pt) = nSr^^y^-^^ {x,pt) = (A9) 



qq-'qq 

J3,(99-99) (^^p^) ^ ^ ^(0) I 3 Tr[^'°'] ^[+] _ 1^P]^[+] _ 1^[+1 _ J ^ (A^O) 

^5^(99-99) (^^p^) OC (^(0) { iwP'w'+' - _ ^ 1^[-] I ^ 



jg,(99--.99-)(^^p^) = 7r5$(?'^-'''(a:,pT) OC (^(0) { _ | ^ (^^2) 



(A13) 
(A14) 
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Sr^"\x,pT) oc {Tr [F(0)Wl+lF(o{i^^i^wHt-iwHt} +F(0)wHF(o{i^^^^ (A15) 

^5r(?p'"(x,pT) oc (Tr [f{0)U^+^F{O{'^ ^'^^™^ U^-^^-^U^-^^}-F{0)U^-^F (A16) 
nSr^^y^^\x,pT) oc (TV [f (0)i^WF(o{iwHt_iZ^,,Wt} + j^(0)wHi.(5){ii,Wt_iZ^^ (A17) 



(5$ 



{x,Pt) 



Tr[U 



■n]t- 



N 



w[+i-iwW}V;(0> 



(AI8) 



gg-'qq 



6V 



{x,pt) 



(99-*99) 
'f 



{x,pt) 



(TV [F(0)wWF(o{i2^^i^wHt_iw[+lt} + j.(o)wHj.(^){i 



(A19) 



(gs->99) 



oc (TV [f (0)t/WF(C){i:^^wHt_i^Ht} _ j^(o)wHf (e){i2M^i,Wt_i^Wt} ] ^ (A20) 



99^99 



5r(")(a^,PT) oc (Tv[F(0)wWF(o{i^5^i^W[-"t-|wWt} +F(0)wHj^(^)|iZ^^ ) (A21) 

n6r'-S'^-''\x,pT) oc (Tr [^F(O^'°'t^WtF(0)Z^'°WW +F(0)Z^Mf (0^'+'^ '^'"^^'°'^ 

+F(0)wHf(o{ ^w[-lt_iZ«wWt }] ) (A22) 



iag^aa) 



{x,Pt) = 



(A23) 
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